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Abstract 



Q^ We present in this paper an approach for computing the homogenized behavior of a 

•^r medium that is a small random perturbation of a periodic reference material. The 

random perturbation we consider is, in a sense made precise in our work, a rare event 
at the microscopic level. It however affects the macroscopic properties of the material, 
and we indeed provide a method to compute the first and second-order corrections. To 
this end, we formally establish an asymptotic expansion of the macroscopic properties. 
Our perturbative approach shares common features with a defect-type theory of solid 
state physics. The computational efficiency of the approach is demonstrated. 
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O 1 Introduction 

K>" Composite materials are increasingly used in industry. For instance, modern aircrafts 

K^ consist, for more than 50%, of composite materials. Generally speaking, composites are 

?-H heterogeneous materials obtained by mixing two phases, a matrix and reinforcements (or 

inclusions). When appropriately designed, these materials outperform traditional mate- 
rials, notably because they combine robustness and lightness. Their use however raises 
new challenges. The behavior of these materials under extreme conditions has to be pre- 
dicted carefully, so as to avoid, in the worst case scenario, separation of the components 
(think of a plane hit by thunder). While it is possible to create an infinity of composites 
starting from the same elementary components, it is out of question to actually construct 
and experimentally test each and every possible combination. Characterizing a priori the 
properties of a given composite material, not yet synthetized or assembled, is therefore 
instrumental. 

A brute force numerical approach, consisting in directly solving the classical boundary 
value problems modelling the behavior of the material, is not practical. The heterogeneities 
indeed often occur at a scale e much finer than the overall typical lengthscale (say, 1) of 
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the material itself. A finite element mesh would, for example, need to be of size less than 
e in order to capture the correct behavior. The number of degrees of freedom would then 
be proportional to e~ (where d denotes the dimension of the ambient physical space) and 
would yield, for e small, a heavy computational cost one cannot necessarily afford. 

The aim of homogenization is to provide a practical alternative to the brute force 
numerical approach. In a nutshell, homogenization consists in replacing a possibly compli- 
cated heterogeneous material with a homogeneous material sharing the same macroscopic 
properties. It allows for eliminating the fine scale, up to an error which is controlled by 
e, the size of this fine scale as compared to the macroscopic size. Homogenization is a 
well-established theory (see |11| for a comprehensive textbook), which, in a simplified pic- 
ture, can be seen as averaging partial differential equations that have highly-oscillating 
coefficients. 

Of course, the structure of the material, and more precisely the way the constituents 
are combined, have a deep influence on the results of the homogenization process. The sim- 
plest possible situation is the periodic situation. At the fine scale, a unit cell is repeated 
in a periodic manner in all directions. Then, in simple cases (say, to be schematic and 
to fix the ideas, linear well-posed equations), the homogenized material is characterized 
only using the solution of simple problems on the unit cell, called the cell problems. The 
role of these cell problems is to encode the information of the micro-scale and convey it to 
the macro-scale. Related cases, such as pseudo-periodic materials, can be treated similarly. 

As Figure [T] shows, real life materials are however not often periodic. In particular be- 
cause of uncertainties and flaws in the industrial process, composites often do not exhibit a 
perfect periodic structure, even though it was the original plan. A suitable way to account 
for this is to use random modelling. Although the mathematical theory for homogeniza- 
tion of random materials under classical assumptions (ergodicity and stationarity) is well 
known, the practice is quite involved. The cell problems are defined over the whole space 
M and not simply on a "unit" cell. The numerical approximation of such problems using 
Monte-Carlo type computations is incredibly costly: the cell problems are truncated on a 
bounded domain, many possible realizations of the materials are considered, averages are 
performed. Consequently, in the context of random modelling, the benefits of homogeniza- 
tion over the direct attack of the original composite material are arguable. 

Our line of thoughts, and the approach we try to advocate here, are based on the fol- 
lowing two-fold observation: classical random homogenization is costly but perhaps, in a 
number of situations, not necessary. A more careful examination of Figure [T] indeed shows 
that albeit not periodic, the material is not totally random. It may probably be fairly 
considered as a perturbation of a periodic material. The homogenized behavior should ex- 
pectedly be close to that of the underlying periodic material, up to a small error depending 
on the amount of randomness present. 

The aim of this paper is to give a practical example of theory following the above 
philosophy. We introduce and study a specific model for such a randomly perturbed 
periodic material, which we also call a weakly random material. More precisely, we are 
interested in the homogenization of the following elliptic problem 

' - dw (^{Aperi^) +brj{^,L0)Cper{^))yu,^ = /(x) in P C M^ 

if = on dV. 




Figure 1: Two-dimensional cut of a composite material used in the aeronautics industry, 
extracted from [IGj and reproduced with permission of the author. It is clear that this 
material is not periodic, yet there is some kind of an underlying periodic arrangement of 
the fibers. 



Here the tensor Ap^r models a reference Z -periodic material which is randomly perturbed 
by the Z'^-periodic tensor Cper, the stochastic perturbation being encoded in the stationary 
ergodic scalar field 6^. In the present work, the law of the random variable brj{x 



IS a 



Bernoulli distribution with parameter 77 (that is, 6^ is equal to 1 with probability rj and 
with probability 1 — r]. Using an asymptotic analysis in terms of r/, we will develop 



an homogenization theory for A,^ 
theory for Aper (f). 
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) = Apf^r (f) + bri (f jw) Cper based on the similar 



,N 



tensors A^' computed on the super cell 
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In short, let us say that the main result of this article is to formally derive an expansion 

A; = A;,, + ijAl+7]^A*+o{ri^), 

where A* and A* are the homogenized tensors associated with Arj and Aper respectively. 
The first-order correction A'l is obtained as the limit, when A^ — )• 00, of a sequence of 

. It is the purpose of Proposition 2 to 
prove the convergence of yl^'". The second-order term A2 is likewise defined as a limit, 
this time up to extraction, of a sequence of tensors A2 when N — t- cxd. The proof of the 
boundedness of the sequence vlg' which implies this convergence up to extraction is not 
given here for it relies on long and technical calculations. Actually, we strongly believe 
that A2 is a convergent sequence and that we can write the expression of the limit. We 
refer the reader to [1] for the details. We also stress that these corrections are achieved 
through purely deterministic computations. 



The above setting is of course one possible setting where we may develop our the- 
ory, but not the only one. More general distributions are studied in |2]. Other forms of 
random perturbations of periodic problems, in the spirit of [7j|, could also be addressed. 
Moreover, we have deliberately considered the simplest possible equation (a scalar, linear 
second order elliptic equation in divergence form) to avoid any unnecessary technicalities 
and fundamental difficulties. Other equations could be considered, although it is not cur- 



rently clear (to us, at least) how general our theory is in this respect. 

With the ideas developped here (and originally introduced and further mentioned in 
[21 |3l m]), we work in the footsteps of many previous contributors who have considered 
perturbative approaches in homogenization. In [15j and [3], a deterministic setting in which 
an asymptotic expansion is assumed on the properties of the material (the latter being not 
necessarily periodic) is studied under the name "small amplitude homogenization". In |14) . 
the case of a Gaussian perturbation with a small variance is addressed from a mechanical 
point of view. Our setting here is particular, because our random perturbation has order 
one in amplitude. It is only in law that the perturbation considered is small. The cor- 
rections obtained are therefore intrinsically different from those obtained in other settings 
(including settings we ourselves consider elsewhere, see IHE]). Also, the present pertur- 
bative theory has unanticipated close connections with some classical defect-type theories 
used in solid state physics. 

We emphasize that, contrary to what is presented in a companion paper for some 
other distributions, the theoretical results we obtain below in the Bernoulli case are only 
formal. We are unfortunately unable to fully justify our manipulations except in the one- 
dimensional case. Nevertheless, we can prove that the terms we obtain as first-order and 
second-order corrections are indeed finite and well defined. Our numerical results, on the 
other hand, show the efficiency of the approach. They somehow constitute a proof of the 
definite validity of our perturbative approach, although we wish to remain cautious. Note 
that due to the prohibitive cost of three-dimensional random homogenization problems and 
the limited computing facilities we have access to, our tests are performed in dimension two. 

This paper is organized as follows. For the sake of consistency and the reader's conve- 
nience, we start by recalling in Section 2 some classical results of periodic and stochastic 
elliptic homogenization. Then we introduce our perturbative model in Section 3, and ex- 
plain how we obtain the first-order and second-order correction by means of an ergodic 
approximation. Our elements of proof are exposed in Section 3. Our two-dimensional 
numerical tests are presented in Section 4. The appendix contains explicit computations 
in the one-dimensional case as well as some useful technical lemmas. 

Throughout this paper, and unless otherwise mentioned, K denotes a constant that 
depends at most on the ambient dimension d, and on the tensors Ap^r and Cper- The 
indices i and j denote indices in [[l,(i]]. 

2 Some classical results of elliptic homogenization 

We recall here some classical well-known results regarding linear elliptic periodic and 
stochastic homogenization. The reader familiar with homogenization theory can easily 
skip this section and directly proceed to Section 3. 

2.1 Periodic homogenization 

Consider A a Z'^-periodic tensor field from R'^ to R'^^'^ , that is 

VA: G Z , A{x + k) = A{x) almost everywhere in x E M . 



We assume that A £ L^{]R. ,]R. ) and A is coercive, which means that there exist A > 
and A > such that 

V^ G R'^, a.e in x e R'^, A|^p < A{x)^ ■ ^ and \A(x)C\ < A|^|. (2.1) 

Consider now a material occupying a bounded domain O C R . The constitutive proper- 
ties of this material are supposed to be periodic, the scale of periodicity being e, and we 
assume that these properties are given by the tensor Ae(x) = ^ (-). 

We consider the following canonical elliptic problem: / € L'^{0) being given, find 
Ue S Hq{0) solution to 



-div(^,Vu,) = /inC', 
Uf = on do. 



f2.2) 



A direct numerical handling of (2.2) using finite elements has a heavy computational 



cost since the scale e of the heterogeneities requires a fine mesh. The aim of homogenization 



is to take the limit e — )• in (2.2) so as to replace the heterogeneous material with a 



homogeneous material. To this end, let us define the periodic cell problems on the unit 
cellQ=[-i,if by: ViG[l,dl, 



- div {A{Vwi + ei)) = in Q, 
Wi Tj — periodic. 



(2.3) 



where et is the i-th. canonical vector of R . Problem (2.3) has a solution unique up to the 



addition of a constant. Note that the number of cell problems is equal to the dimension of 
the space. 



The homogenized tensor A* is then given by: 

y{i,j) G [1, df, A*, = f A{Vwi + ei) ■ ej. 

JQ 



(2.4) 



Using (2.3), it also holds 



A*^ = / A{Vwi + ei)- {Vwj + tj) 
JQ 

Notice that in this periodic setting A* is a constant matrix. 



Finally, let us define the homogenized solution uq as the unique solution in Hq{0) to 

J - div {A*Vuo) = / in O, 
I uo = on do. 



(2.5) 



Solving (|2.3|) and (|2.5|) is much simpler than directly solving (2.2) for the fine scale e has 

(2.6) 



disappeared. It is well-known (see fllj for instance) that 



Ue — > no in L-^iO) 



and 



Ue -Uo 



i=l 



w^A-J^ ^0 \nH\0). 

e/ OXi e-s>0 



(2.7) 



The functions Wi are also called the correctors, since they allow for the strong convergence 



in (2.7). Convergences (2.6) and (2.7) show the relevance of the homogenization process: 
can be replaced by uq or more accurately u^ + e X^j=i Wi (|) ^^(2;), which are easier to 



u. 
compute. 



2.2 Stochastic homogenization 



Throughout the article, {il.,J^,¥) denotes a probabihty space, P the probabihty measure 
and a; € $7 an event. We denote by E(X) the expectation of a random variable X. 

We assume that the group (Z , +) acts on Q and denote by r^, /c E Z , the group action. 
We also assume that this action is measure-preserving, that is, 

and ergodic: 

V^ G T, (VA: G Z'^,^ = TkA) =^ {F{A) = or P(^) = 1). 
We call F G LJ^^iR'^ , L^Q)) stationary if 
VA; G Z , F{x -\- k,uj) = F{x,TkCo) almost everywhere in x G M and uo £ Q. (2.8) 



Notice that the notion of stationarity we use here is discrete: the shifts in (2.8) are 



assumed to be integers. This is related to our wish to connect the random problems con- 
sidered with some underlying periodic problems. Notice also that for a deterministic F, 
stationarity amounts to Z -periodicity. 



Consider a stationary tensor field A{x, u) G L°°{R'^ x i7, M'^^"'), such that (2.1 ) is almost 



surely satisfied by A(-,u}), and a material occupying a bounded domain O C M modeled 
by^(f,w). 

We are interested in solving, for a deterministic function /, 



Mf 



div lAl-,uj] Vu, 
= on do. 



fix) in O, 



(2.9) 



In order to describe the behavior of u^, we again need to define cell problems. Here 
they read (see |llj): 



— div {A{x,CL!){'Vwi + ei)) = in 
Vwi stationary, E ( / Vwi ] = 0. 



(2.10) 



Problem (2.10 ) has a solution unique up to the addition of a (possibly random) constant. 



Then we define the homogenized tensor A* by 

V(i,j) G ll,df,A*, = E(j A{y,u){ei + VyWi{y,uj))-ejdy] . (2.11) 

Notice that A* is deterministic and constant throughout the domain O. The homogenized 



field uq, which gives the asymptotic behavior of u^ (in a sense similar to (2.6) and (2.7)), 
is also deterministic. It is the unique solution in Hq{0) to 

J -div(^*V'Uo) = fi^O, 
I lio = on do. 



The computation of the stochastic cell problems (2.10) is not an easy task since the 
problems are posed in an infinite domain (M ) with a stationarity condition. As we have 



seen in the previous paragraph, when the material is periodic, the cell problems (2.10) 



reduce to the deterministic cell problems (2.3) which are Z -periodic and can thus be 



computed on the unit cell Q. Consequently, when the material under consideration is a 
stochastic perturbation of a reference periodic material, we expect the computation of the 
homogenized tensor to be tractable, up to an approximation. This is our motivation for 
proposing a perturbative approach. 

3 Homogenization of a randomly perturbed periodic material 
3.1 Presentation of the model 



In the stochastic framework (2.9)-(2.10)-(2.11), we now specifically consider the following 
tensor field in W^ x Q: 



Ar^{x,L0) = Aper{x) + br^{x,Uj)Cper{x). 



(3.1) 



Here Aper and Cper are two deterministic Z -periodic tensor fields. Intuitively, Aper 
is the reference material perturbed by Cper- The random character of the perturbation is 
encoded in the stationary ergodic scalar field brj , upon which we assume the expression 

bn{x,uj)= Y^ lQ+k{x)B^{uj), 

where the B^ are independent random variables having Bernoulli distribution with param- 
eter ry, meaning B^ = with probability 1 — rj and B^ = 1 with probability t]. 

It is clear that as r/ — )• the perturbation becomes a rare event. However, the realiza- 
tion of this event modifies the microscopic structure of the material since it replaces, in a 

given cell, Aper with Aper + Cper- 

We additionally assume that there exist < a < /3 such that for all ^ G R and almost 

all X G M*^, 



a|^P < Aperix)^ ■ e, a|eP < {Aper + Cper) (x)^ ' ^, 



(3.2) 



\Aper{xn < m, I {Aper + Cper) {x)i\ < ^V (3-3) 

We can therefore use for every < r? < 1 the stochastic homogenization results recalled 



in Sectional The cell problems associated with (3.1) read, for 1 < i < d, 

' - div {Ar^{Vw1 + ei)) = in M'^, 



VtD^ stationary, E ( / Vw^ j = 0, 



(3.4) 



and the homogenized tensor A* is given by 

A*^ei = ^( I M^^i + ^i)) ' for 1 < i < d. (3.5) 

Throughout the rest of this paper we denote by tf^ the solution to the i-th cell problem 



(2.3) associated with Aper- 



Because of the specific form of Arj, and more precisely because A^ converges strongly 
to Aper in L^{Q X il) as r/ — )• 0, it is easy to see that: 



Lemma 1. When 77 — > 0, A* — > Ap^^. 

Proof. Fix 1 < i < d. We start by proving that Vw^ converges strongly in L^{Q x il) to 
Vw^. Indeed, define r? = w? — w? solution to 



' - div (A^Vr^'') = div (6^Cper (Vu;° + e^)) in 



Vr.' stationary, E / Vr-' = 0. 



(3.6) 



Standard cut-off and ergodicity arguments (see e.g the proof of Proposition 3.1 in j7]) 
show that 

— II rOii iir* fy7o,fi < ^ \ II 

= -VvWCper (Vw° + a) \\l2(Q), 



where a is defined in (3.2), so that Vw^ — )• Vw^ in L'^{Q x Q,) . 



Next, it is straightforward to see that A^j converges strongly to Aper in L^{Q x O). We 
deduce from these two strong convergences that 

A*^ei = ^( Arj{x,uj){Vw'^ + Ci) j ^^ / AperiVw^ + ei) = A*p^^ei. 

This concludes the proof. 

D 

Our goal is now to find an asymptotic expansion for A^ with respect to r] up to the 
second order. 



3.2 An ergodic approximation of the homogenized tensor 

We consider a specific realization u; & Q oi the tensor A^j in the truncated domain 

2 ' 2 J ' 



/tv = [~~2' t] ' "^ith (for simplicity) A^ an odd integer, and solve the following "supercell" 



problem: 

- div [An (x,oj)(Vw'^ ' '"^ + 61)) =0 in Ijy, 

^ ^ (3.7) 

^■n,N,^ (iVZ)'^ - periodic. 

Then an easy adaptation of Theorem 1 of [8], stated in the continuous stationary 
setting, to our discrete stationary setting, shows that when A^ goes to infinity, 

—^ / Ari{x,uj){Vw^' '"^(x) + ej)(ix converges to A*ei almost surely in a) G 0. (3.8) 
J In 

Since —7 / An{x,oj){'S/w^' ''^(x) + ei)(ix is the tensor obtained by periodic homogeniza- 

•^ In 
tion of the tensor A^(a;,a}) in the supercell /jv, it is also well-known (see |llj ) that the 

following bounds hold for all {i,j) £ [[l,dp: 



-^ ( / Aj^^{x,u))dx] et-ej < --^ A^(x,tIJ)(Vu;f ^''^(x) + ei) ■ ejdx 



As a result, for all N in 2N + 1, for all < ry < 1 and almost all w in fi, 



r),N,Cj I 



—^ A^(x,(Ij)(V<' ' {x) + ei)-ejdx 
J In 



</3, 



(3.9) 



where /3 is defined in (3.3). We then deduce from (3.8), (3.9) and the Lebegue dominated 



convergence theorem that 



\d^(l A^(x,a;)(Vu;f''^'"(x) + ei)')dx. (3.10) 



Ml < i < d, A*ei = lim ,^ , 

Remark 1. A similar result holds for homogeneous Di rich let and Neumann boundary 
conditions instead of periodic conditions in the definition (3.7) of w^' '^ (see lE^ for more 
details). 

Using now the fact that brj has a Bernoulli distribution in each cell of Z , it is a simple 



matter to count the events and to make (3.10) more precise. We first define the set 

d 



Tn = { 



k£Z'',Q + kcl 



N 



N-1 N-1 



(3.11) 



The cardinal of T/v is of course A^ , and M {Q + k} = I^. 

keTN 
We then have the following possible values for j4.^: 

• Ari{x,Cd) = Aper with probability (1 — r/) . 

In this case w^' '^ = w^ solves the usual periodic cell problem: 



div [AperiVw^ + a)) =0 in Q, 



uij Z — periodic. 



JV^-l 



In this case w^'"''^ = tf^''^'" solves the following problem, which we call here a "one 



Ar^{x,oj) = Aper + l{Q+fc}Cper for k £ Tn, with probability 7/(1 — rj) 

ri.N.u l,k,N 

defect" supercell problem: 

l,fc,7V 



div ( [Aper + l{Q+fe}Cper) (Vw- ' ' + ej) j = in I 



'Af, 



w. 



l,k,N 



{NZ)'^ - periodic. 



(3.12) 



Arj{x,Cd) = Aper + '^{Q+i}u{Q+m}Cper for (l,m) £Tn, I y^rn, with probability 

In this case w^' '^ = w-' '"*' solves the following problem, which we call here a "two 
defects" supercell problem: 



div ( {Aper + l{Q+i}U{Q+m}Cper) (Vw- ' ''"' + Cj) ) =0 in I^, 



2,l,m,N 



{NZf - periodic. 



(3.13) 



All the other possible values for ^^, which are of probability less than rf and which 
we will not use in this article, can be obtained using similar computations. 



An instance of a setting with zero, one and two defects is shown in Figure |2] in the 
two-dimensional case of a material Ap^r consisting of a lattice of inclusions. 



• ••«••••# •♦•♦ •«•• •«•« •♦•• 



Figure 2: From left to right: zero defect, one defect and two defects. 
Let us define A^ = Aper + l{Q+fc}Cper and vlg'™ = ^per + l{Q+i}u{Q+m}Cper- Then 



( |3.10[ ) reads 



^^^ / A^ (Vwr +e,) + 



It . clear, by (iVZ)^-pen„d.ity, that / ^J(V».'*- + e.) doe. not depend on the 

position k G T/v of the defect. Likewise, / A2^{'S/w-' '™' + Cj) only depends on the 

Jin 
vector m — I. Thus we can rewrite 

A*e. = hm f (1 - 7?)^'A;,,e. + r/(l - r/)^'-i / A?(V^1'°'^ + e.) 

+ Y. ""'^^-f'' [ ^f(Vt.^°'^'^ + e.) + •••). (3.14) 

fee' 

This is of the form 



k£TN,k^O ■'^^ 



:,Af 



^* = lim y r/PA!'- 

p=0 (3.15) 






where the remainder Oiv(^^) depends on A^. 



Explicitly expanding the polynomials in rj up to the second-order in (3.14), we obtain: 



^o' ~ Ap^^, (3.16) 

A*'^e,= / A\{Vw]^''-'' + e,)- j Vr(V«;0 + e,), (3.17) 

J In JIn 



10 



^2 ^* 



2 Z^ [ ^2 i^w-' ' ' +ei)-2 A^iVw^' ' + ei) 



feeTiv.fc^o 



(3.18) 



as the first three coefficients in (3.151. 



Remark 2. The structure of Ap for p S N is obviously related to that of the polynomial 
(l-x)P. 



Our approach consists in formally exchanging the limits N ^ oo and r/ —t- in (3.15). 
In the next section, we show that ^4^' is a converging sequence when N — )■ oo. The case 
of A2' , which is shown to be a bounded sequence and thus to converge up to extraction, 
is discussed in Section [3.41 



.*,N 



J2 A*,N-s 



We are not able to prove, though, that A* — lim {A* — r]A*^ — r/ vlg' ) = o{r] ) 

with a remainder term o{rf') independent of N . 

Remark 3. The expression of A*^ (and likewise A*2 ) is reminiscent of standard ex- 
pressions in solid state theory: each of the two integrals in the definition (3.11) of A*^ 



scales as the volume N of the domain Ij\f, and a priori needs to be renormalized in order 
to give a finite limit. The difference however has a finite limit without renormalization. In 
solid state physics, it is common to substract a jellium, that is, a uniform background, and 
proceed similarly. 



3.3 Convergence of the first-order term Al 



,N 



N 



We study here the convergence, as N goes to infinity, of A^' defined by (3.17), and prove: 



Proposition 2. yl^' converges to a finite limit A* in M when N — t- 00. 



*,N 



Proof We fix {i,j) in [1, dj and study the convergence of A-^^" a ■ e 



Let us define the adjoint problems to the cell problems (2.3): 



- div (A^^^(Vu;° + Cj)) =0 in Q, 
Wj TL — periodic, 



(3.19) 



where we have denoted by ALg^j. the transposed matrix of A^er- Then using (3.12) and the 



definition of Al[^ we have 

/ ^?(Vu;^°'^ + e,).e, = / A?(V^1'°'^ + e.) • (e, + Vu;°) 

= f Aper{Vwl''''' + ei)-{ej + Vw']) 
J In 

+ [ CperiVwl''^''' + ei) ■ {ej + Vu;°). 
JQ 



Next, using (3.19), we note that 



f AperiVwl''''' + ei)-iej + Vw'^) = f (V«;^°'^ + e,) • Ajer(ei + V 
Jlisr Jiff 



»«) 



/ 

JIn 



ei-^L(ei+Vu;°) 



^ [Aper) ^j ' Sj) 



11 



and applying (2.4) to the periodic tensor Ai^j. and noticing that {A±^^)* = (A*), we 



obtain 



[ A^^iVw]^'"'^ + et) ■ e,- = N^A^^^Ci • e,- + /" Cper(Vu;,^'°''^ + a) ■ (e,- + Vu;°). (3.20) 
Jin Jq 

Since, by definition, 

J In 'I In 

= / Ai{\/Wi' +ei 

JIn 



'In 

N*^ A* P- 



we deduce from (3.20) that 






Jq 



We now define 



l,0,Af l,0,Af 



whicli solves 



- div (^?Vg^°'^j = div(lQCper(Vu;0 + e,)) in In, 
gi'O'^ (ATZ)'^ - periodic. 



(3.21) 



(3.22) 



(3.23) 



We deduce from Lemma k)^ of the appendix, applied to (3.23), that Vg^' ' converges in 

Lf^^{]R.'^), when A^ — )• +oo, to Vq-' '°°, where q-' '°° is a Lf^^{]R.'^) function solving 



- div (a;v(?^°'°° 

^1,0,00 ^ ^2(j^d)_ 



div {1 qCper C^ Wi +ei)) in 



(3.24) 



Defining w. 



1,0,00 



I 1,0,00 



It follows from (3.21) that Al 



,*,N 



N^+cxi 



lfi,N 
i 

A\ with A\ defined by 



w^ + q^' '°°, it is clear that Vit^j ' ' converges in L {Q) to Vif, 



1,0,00 



V(i,i)G[l,df, A\ei-e, 



1,0,00 



Cper(Vz/;,^'"'°° + Ci) • (e,- + Vu;^^). (3.25) 



D 



The computation of A\ requires to solve (3.24) which is defined in R , but, in sharp 



contrast to the stochastic cell problems (3.4), is deterministic and has a right-hand side 
with compact support in M . In practice, problem (3.24) is truncated on I^- The following 
result gives insight on the truncation error. 

Lemma 3. Assume that d > 3 and that the unit cell Q contains an inclusion D, the 
boundary of which has regularity C^''^ for some < fi < 1, and such that dist{D,dQ) > 0. 
Assume also that Ap^r is Holder continuous in D and in Q\D. Then there exists a tensor 
B^' , computed on Ijv, and a constant K independent of N such that 
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Proof. Step 1. 

Fix {i,j) in |{l,(ip. We first define the adjoint problem for (3.23), namely 



div ( {A'ifVq]''''' 



dw{lQC^^^.{Vw'> + ej)) in In 



~1,0.N 



~i,U,iV ^j^^y _ pg^iQ^ip_ 



(3.26) 



^j 



Applying Lemma 6 to (3.26 



It solves the adjoint problem o 



-1,0,00 f ~lfl,N 



, we also introduce the limit q-' ' of q,-' ' when N ^ oo. 
(I3^. 



Then, using (3.26), we obtain 



JQ JIn 

f ^Ov7 l,0,Af ^~1,0,N 
J In 



1,0, N 



Consequently, (3.21) and the definition (3.22) of g^' ' yield 

Al'^'e.-ej = fCper{Vwf + e,)-{e, + Vw^)-[A^,Vql^''''-Vq]^^'''. (3.27) 
jQ JIn 

We knowfrom Lemmap applied to (3.23) and (3.26) that 1/j^Vg^' ' andl/j^Vq"' ' converge 



roToo 



strongly in L (M ) to Vg^ ' ' and Vg ■' ' respectively, when N — )• (X). Passing to the limit 
in (3.27) then gives 

1^* ■ Cj ^ 



Cper(Vw^ + ei)-{ej + Vw°j 



AO-r-7 1,0,00 r-7~li0,OO 



We now define v^ ' ' andv' ' solutions to (3.23) and (3.26) with homogeneous Dirich- 
let (instead of periodic) boundary conditions on the boundary OIn of I^, and the ten- 
sor n^ by 



-t>l Cj • €j 



[ CperiVw^ + ei) ■ {ej + Vw^A - [ A?Vt;,^'°'^ • Vv 
JQ JIn 



3 



(3.28) 



The proof of Proposition ^ is easily adapted to show that B\' converges to A\ as N 
goes to infinity. 

Step 2. 



We consider 



-D]^ A-^ I €-1 • €.j 



/|0v7 1.0,00 v^~l,0,OO / 

Jl 



.Ov7 l,0,Af v7~1.0,Af 



In 



-,*,N 



and expand the difference B^' — A^ as follows: 



Bl'^-Al]e..e, 



AKJT7 -L,U,CXJ VT~-L- 

A^Vq^' ' -Vq- 



I I I /lOvT 1,0,00 r-T~li0,OO 



f aOv7 1:0-^ V7~1,0,A'A 



(3.29) 
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We now show that the two terms in the right-hand side of (3.29) converge to as A^ 
when N — )• +00. 

We first note that the results of Lemma [s] of the appendix, stated for a Z'^-periodic 
matrix, can be readily extended to address A^ since A^ is equal to Aper in M. \Q. 

We deduce from Lemma 7 applied to (3.24) that g^' '°° is defined uniquely up to an 
additive constant. Moreover, Aper being piecewise Holder continuous, we deduce from 
Lemma Is] that there exists a unique solution to (3.24) which converges to zero at infinity. 

Since we only use Vq- ' '°° in A^ , we can thus assume without loss of generality that 
q^' '°° converges to zero at infinity. Likewise, we assume that q-' '°° converges to zero at 
infinity. 

We then deduce from Lemma [8] that there exists a constant K independent of A^ such 
that for Ixl > 1, 



k-'°'°°(x)| < K\x\^-'', |^j'°'°°(x)| < i^|x|i-^ 

|Vg,^'°'°°(x)| < K\x\-'^, |VqJ'°'°°(x)| < K\x\-'^, 

IvF'^'ixM < K\x\^-^ |t}l'°'^(x)| < K\x\'-'', 



J 



lfl,N , 



-lfi,N, 



Using (3.31), we have II Vg^' '°°\\i2 



and 



Vvl'""''' {x)\ < Klxl'"^, \Vvp''{x)\ < K\x 
n) 



I /lOv7 1,0,00 v7~l 






<KiV-^/2and||Vgj'°'°°||i2 

l,0,OO|| ||VT~1,0,00|| 



lL2{]Rd\7iv)ll^^i 



lL2 



(3.30) 
(3.31) 
(3.32) 
(3.33) 



(3.34) 



where /3 is defined in (3.3). 



We now address the second term of the right-hand side of (3.29) and write 



.Ov7 l,0,Af V7~l,0-Af 



/lOvT 1,0,00 r-T~l,0,OO 



f /lO/v7 l,0,Af V7 l,0,ooN v7~l,0,Af , f hOvt 1,0,oo f^~lfl,N ^7 ~ 

J I AT J IfJ 



1,0, oo\ 

i > 



Since div 



l,0,Af 



1,0,00- 



vr- -Vql^'^n) = div((yl'^)'^(V{)j'°'^-Vg]'°'°°) 



~' ' ' = on 9/7V, we have, using integration by parts. 



^j 



V7 l,0,oox ~l,0,Af 



f /lO/v7 l,0,iV V7 l,0,ooN ^~lfi,N , I .Ov7 1,0,00 /v7~1.0,Af V7~l 

/ ^i(VVi'' -Vg^'- )-V7;/' +/ AiVg^' • -(Vf^' -Vq^ 

Jiff JljSi 

= [ A?(V.1'°'^ 

JdiM 

+ / (^?)^(V^]'°'^ - v,~ 

JoIn 

f I /lO\T/v7~l,0,Af v7~l,0,ooN 
= / (^1) (Vv/ -Vg/ )-vq^^ 
J SIn 



in In, and 



1,0,OON 



-1,0, oo\ 1,0,00 



1,0,00 
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where v is the unit outward normal vector to dl, 



N- 



The estimates (3.30) and (3.33) imply 



I 1,0,00 II 

m IIl°°(9/jv) 



<KN'~'', \\{A^,fiW'''^ 



while the measure of the boundary OIn scales as N . Hence 






Vg~J'°'°°) • vU^^ai^) < KN"", 



Join 



~l,0,OO\ 1,0, CXD 



< KN- 



and then 



.Ov7 l,0,Af v7~l:0,Af 



/lOvT 1,0.00 r-T~l,0,OO 



< KN- 



We conclude by substituting (3.34) and (3.35) into (3.29). 



(3.35) 



D 



Remark 4. We assume d > 3 and piecewise Holder regularity on Aper, and use Dirichlet 
boundary conditions in Lemma^ because our proof relies on Lemma^ Note however that 
the numerical experiments of Section^^show, in dimension d = 2, that we again obtain 
the rate N~ in the convergence of A*^ to A\ for two different Ap^.^, one being piecewise 
Holder continuous in the sense of Lemma^and the other not, and with periodic boundary 
conditions. Moreover, the explicit computations of Proposition [^ show that in dimension 
, and without any assumption of regularity on A^, the rate of convergence of A'l' to 



one 



Al is iV-i. 



,*,Af 



3.4 Second-order term 

For completeness, we state here the result regarding A2 proved in [T]: 

Proposition 4. The sequence A^^ defined by (3.18) is bounded in W^^'^ and therefore 
converges up to extraction. 

We strongly believe that A*2 is actually a convergent sequence, as shown by our 
numerical tests thereafter. In fact, we even believe (see [1]) that we can write the expression 



of the limit. Note also that the explicit computations of Section 5.1 prove the convergence 



of A*2 in dimension one. 



4 Numerical experiments 

Our purpose in this section is to assess the approximation of A* by the second-order 
expansion A*^ + tjA^^ + rf^ ■ The limited computational facilities we have access to 
impose that we restrict ourselves to the two-dimensional case. We first explain our general 
methodology and then make precise the specific settings. 

4.1 Methodology 

We will consider two commonly used composite materials as periodic reference materi- 
als Aper- The first material consists of a constant background reinforced by a periodic 
lattice of circular inclusions, that is 

Aper{xi,X2) = 20 X Id + 100 ^ T^B{k,0.3)ixi,X2) X Id, 

k&I? 
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where B{k, 0.3) is the ball of center k and radius 0.3. The second material is a laminate 
for which 

Aper{xi,X2) = 20 X Id + 100 ^ lKa;i<i+l(xi, X2) X Id. 

l&Z 

In the case of material 1, the role of the perturbation is, loosely speaking, to randomly 
eliminate some fibers: 

Cper{xi,X2) = -100 ^ lB{k,0.3)(.Xl, X2) X Id. 

keZ'^ 

In the case of material 2, the perturbation consists in a random modification of the lami- 
nation direction: 

Cper{xi,X2) = -100'^li<xi<i+i{xi, X2) X Id + I00'^li<x2<i+i{xi, X2) X Id. 

lez lez 

In both cases, we have chosen the coefficients 20 and 100 in order to have a high con- 
trast between Ap^r and Ap^r + Cper, and thus for the perturbation to be significant. There 
is of course nothing specific in the actual value of these coefficients. 

These two materials are shown in Figure [3j 



Figure 3: Left: a periodic lattice of circular inclusions. Right: a one-dimensional laminate. 

Our goal is to compare A* with its approximation A* + rjA^' + ff'^ for each of 
these two particular settings. A major computational difficulty is the computation of the 



"exact" matrix A* given by formula (3.5). It ideally requires to solve the stochastic cell 
problems (3.4) on M . To this end we ffist use ergodicity and formulae (3.7) and (3.10), 



and actually compute, for a given realization w and a domain I^ which is here equal to 



*,N , 



[0, A^] for convenience, Arl (w) defined by 



A*^'\u)ei 



1 



Ari{x,uj){Vivl' ''^(x) -I- ei)dx. 



In a second step, we take averages over the realizations co. 



(4.1) 



For each oj, we use the finite element software FreeFem++ (available at www.freefem.org) 



to solve the boundary value problems (3.7) and compute the integrals (4.1 ). We work with 
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standard PI finite elements on a triangular mesh such that there are 10 degrees of freedom 
on each edge of the unit cell Q. 



t,N 



■■,N , 



We define an approximate value j4JJ'^^ as the average of A^' (cj) over 40 realizations u. 
Our numerical experiments indeed show that the number 40 is sufficiently large for the 
convergence of the Monte-Carlo computation. We then let N grow from 5 to 80 by incre- 



ments of 5. We observe that A^' stabilizes at a fixed value around A'^ 
Ari' as the reference value for A* in our subsequent tests. 



80 and thus take 



The next step is to compute the zero-order term A*^^ and the first-order and second- 



Using the same mesh and finite elements 



order deterministic corrections A*^' and ^2' ■ 

as for our reference computation above, we compute j4Lj. using (2.3) and (2.4), and for 

each A^ we compute A^ and A2 using (3.17) and (3.18). We again let N grow from 



t:,N 



:,Af 



5 to 80 by increments of 5 for A-^ . The computation of A2 being significantly more 
expensive (note that in (3.18) there is not only an integral over /^v but also a sum over 
the iV^ cells) we have to limit ourselves to A^ = 25 and approximate the value for A^ larger 
than 25 by the value obtained for A^ = 25. 



Before presenting our results, we wish to discuss our expectations. Note that there are 
three distinct sources of error: 

• the finite elements discretization error; 



the truncation error due to the replacement of M with /tv, in the computation of 
the stochastic cell problems (3.4) that are replaced with (3.7), as well as in the 



computation of the integrals (4.1 ); 



• the stochastic error arising from the approximation of the expectation value (3.10) 
by an empirical mean. 

The discretization error originates from the fact that, in practice, we only have access 
to the finite element approximations of all the functions manipulated here (such as w?, 
w^' ''^,...). Although we have not proved it in the specific context of our work, we believe, 
because it is shown in a similar weakly random setting (see [U]), that all the convergences 
stated here in the infinite-dimensional setting still hold true for the finite-dimensional ap- 
proximations of the objects. Our numerical results indeed confirm it is the case. In order 
to eliminate the discretization error from the picture, our practical approach consists in 
adopting the same finite element space for all approximations of the cell and supercell 
problems, independently of A^. 



The truncation error is a different issue. For the "exact" computation of A* (we mean 
not using the second-order expansion (3.15), but (4.1)), we use an empirical mean and a 
truncation. We know from |8], for a continuous notion of stationarity analogous to the 



discrete notion (2.8) we use here, and under mixing conditions which are satisfied in our 



setting, that the convergence of the truncated approximation to the ideal value holds at 
a rate A'^"'^ with At a non explicit function of the dimension, the mixing exponent and the 
coercivity constant of the material. On the other hand, in the second-order expansion 
(3.15), the zero-order term A*^ is of course free of any truncation error. All that we know 



.*,N 



for the approximation A^ defined by (|3.17) to the first-order correction A\, is stated 



in Lemma [3] in dimension d > 3, under Holder regularity assumptions on Ap^r^ and with 
Dirichlet boundary conditions replacing periodic ones. One of the aims of our experiments 
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is therefore to draw some numerical conclusions on the convergence of this term when these 
assumptions are not satisfied. Note that the matrices involved in our test materials are 
clearly discontinuous functions of x. The matrix corresponding to material 1 is piecewise 
Holder continuous in the sense of Lemma [31 while the matrix corresponding to material 2 
is not. As for the second-order approximation vlg' , we have no insight on the truncation 
error and we also wish to study its convergence from a numerical point of view. 

Finally, we have a practical approach to the stochastic error: besides the empirical 
mean, we provide, for each A^, the minimum and the maximum values of A^' (cj) achieved 
over the 40 computations. 

We now would like to emphasize that the purpose of our numerical tests is not to prove 
that 

a; = 4,, + r?^t'^ + ifA*^"^ + o{r^') 

for a remainder term o{rf) that is independent of N ^ of the number of realizations and 
of the size of the mesh. Establishing experimentally that such an asymptotic holds is too 
demanding a task. It would indeed require letting r] go to 0, which in turn, since we have 
to observe at least one (and in fact many) event per domain considered, would necessi- 
tate a supercell of size N extremely large. We cannot afford such a computational workload. 

Using our numerical tests, we only hope here to demonstrate, and we indeed do so, that 
the second-order expansion is an approximation to A* sufficiently good for all practical 

purposes, and in particular for r/ not too small ! We will observe that A*^ is not only 
bounded as stated in Proposition Wl but, as A^ goes to infinity, converges to a limit A2, 
and that both A*^ and A*2 converge to their respective limits faster than A^' to A* 
(which is intuitively expected since the former quantities are deterministic and contain less 
information). We will also observe that A*^. -\- r]A*^ is significantly closer to A* than 
A*, thereby motivating the expansion. The inclusion of the second-order term further 
improves the situation. 

4.2 Results 

In order to give an idea on how the perturbation affects the materials considered, we first 
show some typical realizations in Figure [4] and Figure [5] Our results are presented in 



Section 4.2.1|and Section 4.2.2 below. Since these results are qualitatively similar for the 



two materials, we comment on the results altogether in Section 4.2.3 
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Figure 4: Two instances of material 1 with t] = 0.1 (left) and ij = 0.4 (right). 
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Figure 5: Two instances of material 2 with r/ = 0.1 (left) and r/ = 0.5 (right). 

To present our numerical results, we choose the first diagonal entry (1,1) of all the 
matrices considered. Other coefficients in the matrices behave qualitatively similarly. As 
mentioned in the previous Section, we illustrate a practical interval of confidence for our 
Monte-Carlo computation of A* by showing, for each A^, the minimum and maximum val- 
ues of A^' (oj) achieved over the 40 realizations oj. 

We will use the following legend in the graphs: 

• periodic: gives the value of the periodic homogenized tensor A*; 

• first-order: gives the value of A* -\- rjA^' ; 

• second-order, gives the value of A* -|- r]A*^' + rj^A^ ; 



*,N 



• stochastic mean, minima and maxima: respectively give the values of A^q' and the 
extrema obtained in the computation of the empirical mean. 

Finally, the results are given for some specific values of r/ (not necessarily the same for 
both materials) which serve the purpose of testing our approach in a diversity of situations, 
from a "small" to a "not so small" perturbation. 
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4.2.1 Results for material 1 

We show the results for ry = 0.1, r] = 0.4 and rj = 0.5 (Figures pi ^ and ^ respectively) . 
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Figure 6: Results for material 1 and r/ = 0.1. Above: complete results. Below: close-up 



on A^' and the first and second-order corrections. 
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Figure 7: Results for material 1 and r] = 0.4. Above: complete results. Below: close-up 
on Ajj' and the first and second-order corrections. 
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Figure 8: Results for material 1 and r] = 0.5. Above: complete results. Below: close-up 
on Ay and the first and second-order corrections. 
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4.2.2 Results for material 2 



We now show for material 2 tlie results for r/ = 0.1, rj = 0.3 and rj = 0.4 (Figures^ 10 and 
11 respectively). 
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Figure 9: Results for material 2 and 7] = 0.1. Above: complete results. Below: zoom on 
Ari' and the first and second-order corrections. 
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Laminate- eta=0.3 
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Figure 10: Results for material 2 and rj = 0.3. Above: complete results. Below: zoom on 



A^' and the first and second-order corrections. 
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Laminate- eta=0.4 
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Figure 11: Results for material 2 and rj = 0.4. Above: complete results. Below: zoom on 



A^' and the first and second-order corrections. 
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4.2.3 Comments 

Notice on the results for both materials (it is especially clear on the close-ups) that the 
first and second-order corrections A^' and Ag' converge very fast in function of N, and 
in particular, as expected, much faster than the stochastic computation. Convergence of 
these deterministic computations is actually typically reached for A^ = 10. 

Then, for all values of rj, it is clear that the first-order correction enables to get sub- 
stantially closer to A* The interest of the second-order term is also obvious as r/ gets 
larger, and we stress that the results are still excellent for r] as large as 0.5, so that our 
approach is robust. 

It is interesting to get some insight on the rate of convergence of the first-order correc- 
tion, and to see whether the theoretical results of Lemma [3] still hold beyond the somewhat 
restrictive assumptions set in this lemma (d > 3, piecewise Holder regularity on Aper and 
Dirichlet boundary conditions on OIn). Recall that d is equal to 2 in our tests, and that 
A'l' is computed with periodic boundary conditions on the supercell /at. Moreover, while 
the lattice of inclusions is piecewise Holder continuous in the sense of Lemma [S] (meaning 
that there is an inclusion stricly contained in the unit cell Q and that the matrix Aper is 
Holder continuous in each phase), the laminate is not. 

We thus plot, for A^ going from 1 to 20 and for both materials, log(| {A^' — A\)ei • ei|) 
in function of log(A^). We recall that A\ is numerically given by ^^' . For both materials 



the 20 points are arranged in a straight line (Figures 12 and 13). This leads us to perform 



a linear regression in order to obtain the slope of the lines. As regards material 1, we find 
a slope of —2.05 and a coefficient of correlation R = 0.99. For material 2, the slope is —1.9 
with a coefficient of correlation equal to 0.95. The rate of convergence for both materials is 
then approximately 0{N^'^) with d = 2, which seems to indicate that the result of Lemma 
|3] still holds true in these circumstances. 
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Figure 12: Rate of convergence of the first-order correction for material 1. 



26 









La 


min 


ate - Rate of converg 


ence 


of the first-order correction 


0.5- 
















0.0- 




X 












0.5- 










X 






1.0- 










X 






1.5- 










X 


















X 




2.0- 














X 


2.5- 










slope =-1 .905492 




X 
X 


3.0- 
3.5- 
4.0- 










R =0.9522175 




X 
X 





















5 








1 .0 1 .5 




2.0 2.5 3. 



log(N) 



Figure 13: Rate of convergence of the first-order correction for material 2. 



5 Appendix 



The purpose of this appendix is two-fold. In Section 5.1 we prove that the approach exposed 
in Section |3] which rehes on formal considerations for general dimensions, is rigorous in 
dimension d = 1. In Section 5.2, we prove for convenience of the reader some technical 
results used in Section |3l 



5.1 One-dimensional computations 

Although we are aware that homogenization theory is very specific in dimension 1 , and can 
be somehow misleading by its simplicity, it is still important to check that our approach 
is rigorously founded in this setting. This is the aim of this section. 

To stress that we work in dimension one, we use lower-case letters Cper and Cper instead 
of Aper and Cper, respectively, as well as for all the tensors manipulated. 



We recall that in dimension one, a* and a* are given by the explicit expressions 

-1 / .1 \ -1 



'"per 



1 



E 



i ar 



1 



Ji ^Tjer I Urj^' 



per 



This enables us to prove the following elementary result which shows that our approach is 
correct in dimension one: 

Proposition 5. In dimension d = 1, it holds 



at = a*p^^ + r]dl + rfdl + 0{rf) 

*,N J *,N 



"r) 



where a| and a2 are the limits as A^ —)• oo of a^' and Og' defined generally by (3.11) and 



(3.18) respectively. 
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Proof. We compute 



(«;) 



*\-i 



1 — r/ 



+ 



1 \ a. 

1 

2 1 



per 



Ojper I Cper 



1 a. 



+ r] 



per 



i XOiper ~r Cper 



-'per 



— yO'per) I J^ VO-per 

This yields the expansion 

* * I ^ * \2 

'^rj ~ '^per + V\^per) 



^per 



G'per\P'per ~r ^per) 



"per 



1 OperyO'per + Cpgrj 



+ ^'(a;er)^ 



C, 



•per 



— i '^per V^per ~r C 



+^'(«;er)' 



'^per + ^v'^per J 



L-per 



1 CLpgr yO'per ~t~ ^perj 



1-r/a! 



~'per 



1 ttperyOper \ Cper , 



per 



+ rf{aler? 



-vHalerf 



jiper y^per 
C-per 



per ) 



-1 



1 OperyOper y ^perj 



L^per 



— ^ Ojper \Ciper ~r C 

1 

2 Cper 

1 OiperyOiper ~T Cper) 



per J 



s- 



It follows from (3.9) and (3.10) that the function r/ — t- a* is bounded on [0, 1]. Therefore 



a* =a* + 7]{al 



a* )'^ 

"rj ^per ' 'iK^per) 



^per 



— i O-peryCLper -\- Cper) 



(5.1) 



2/ * \3 



+ r] {a* ) 



^per 



Ojper\^per ~r C 



+ 0(r/3). 



-per^ 



We now devote the rest of the proof to verifying that the coefficients of rj and r/^ in 



(5.1) are indeed obtained as the limit as A^ — )• oo of o*' and Og' generally defined by 



(|3.17) and (3.18) respectively, in this particular one-dimensional setting. 

The one-defect supercell solution w^'^'^ generally defined by (3.12) satisfies here 



whO,N ^ _ pgriodic. 



iV N 



We easily compute 



<,!(A„.o,«,„ 



N 



N 



N_ Cjjer "I" lr_l UCper 



jiper 



Nialer)-' - l\ 



-1 



■^per 



-1 



OperyOper ~r Cper) 



^ I 
^per "•" 



y'^per) 

N 



L,per 



1 O^peryOjper \ Cper) 



+ 



a* )^ 

^perJ 



'-per 



- \ OjperyOper ~r C 



-per J 



+ o{N- 
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*,N 



Thus a{ defined generally by (3.17) takes here the form 

/V 1 

d 



:,Af 



N dx 

2 



l,0,Af 



+ 1 - Na 



= (a* 1^ 

per V^per J 



^per 



and 



*,A'' -* / * \2 

1 0>per\^per ~r Cper) 



Likewise, for fe G [- 

0,kf d _2,0,k,N 



JV-1 Af-l ji 
2 ' 2 II' 



1 CLper\(^per ~t~ C-per, 



-'per 



+ o{l), 



(5.2) 



"2 ^^^ 



+ 1 



Af 



N_ 
■ 2 



'-f "Pe'- + l[-|,|]U[A:-|,A:+|]^Pe^', 



per 
* \-l 



N N{a*) 



uper 



={< r) \per) I ^ Cper 



_ 1 dperx^per ~r Cpej- 

+ 4 ^ 'P^^' 



fa* )3 "1 



Ar2 



uper 



OjperyO'per ~r CpgrJ 



+ o(iV" 



which is independent of k (and so of the distance between the two defects). Hence, 



a*2 defined generally by (3.18) writes here 



^2 ~ K'^per) 




Cper 



CiperyO'per ~r Cp^r) 



+ o{l), 



and 




(5.3) 



Using (5.1), (5.2) and (5.3), we verify that 

al = a* + r]al + r]'^al + 0{rf) 



U 



Remark 5. The fact that the distance between two defects does not play a role in the 
computation o/og' is of course specific to the one- dimensional setting. As we have seen, 
this is not true in higher dimensions where the geometry comes into play. 



5.2 Some technical lemmas 

The second part of this appendix is different in nature. We prove here three technical lem- 
mas that are useful for our proofs in Section |3] These results, or related ones, are probably 
well known and part of the mathematical literature. We prove them here under specific 
assumptions for the convenience of the reader and for consistency. We acknowledge several 
instructive discussions with Xavier Blanc on the content of this section. 



We recall that Q 



.1 11 

2' 2J 



and Ij\f 



N N]d 
2 ' 2 J ■ 
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Lemma 6. Consider f £ L'^{Q), and a tensor field A from, M to M such that there 
exist A > and A > such that 

Ve G M'^, a.e in x G R'^, Aj^j^ < A{x)C • C and \A{x)^\ < A|^|. 



Consider a^ solution to 



- div {AVq^) = div(lQ/) in In, 
q^ {NZy - periodic. 



(5.4) 



Then Ij^Vq converges in L {M. ), when N goes to infinity, to Vq°° , where q°° is a Lf^^{l 
function solving 



- div {AVq°°) = div(lQ/) in : 



(5.5) 



Proof. We first obtain a bound on ||Vg ||j;^2(/^\ and then, by compactness, extract a limit 
of tliis sequence. 



Multiplying the first line of (3.23) by q and integrating by parts yields 



N T7„N 



AVq'' ■ Vq 



f-^q 



N 



Q 



from which we deduce 



^q''\\LHi.)<^\\f\\LHQy 



A' 



(5.6) 



(5.7) 



Consider now a bounded domain V C M . For N sufficiently large, we have V C In 
and so 



1 



\\'^Q''\\LHV)<j\\f\\LHQ)- 
Thus Vq is bounded in L'^{T>) for every bounded subset V C M . 

Using diagonal extraction and the weak compactness of Lf^^iM. ), we can classically 
find a subsequence of Vq^ such that, without changing the notation for simplicity. 



Vq'^ -^ h weakly in Lf^M)- 



(5.8) 



We deduce from (5.7) and (5.8) that for every bounded subset D (ZW^, 

Ml^v) < \\\f\\LHQy 

This implies that the vector h is in L^(]R'^). 



We also deduce from (5.8) that for ail {i,j) e ll,df, p^ = §^. This implies that h is 



dxi dxj 
the gradient of a function we call q°°. Since h G L^(M'^), Vg°^ = h is in L^(M'^) and q°° in 



Finally, (5.8) yields that Vq converges to Vq°° in 2?'(M ). We can then pass to the 



limit A^ — 7- oo in the first line of (5.4) and obtain 



-div {AVq°°) = div(lQ/) in M"' 
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in the sense of distributions. 



where q°° is in Lf^^{Mr) and solves 



We have proved that Vq converges up to extraction and weakly in Lf^^{M. ) to \7q°° , 



(5.9) 



- div {AVq°°) = div(lQ/) in M'^, 



We deduce from Lemma ^ thereafter that (5.9) has a solution unique up to an additive 
constant, so that Vq°° is uniquely defined. A classical compactness argument then yields 
that the whole sequence Vq converges weakly to Vg°° in Lf^^{M. ). 



It is clear from what precedes that 

1/ivVg^ ^ Vq°° weakly in L'^iR'^). 



(5.10) 



We now prove that the sequence Ijj^Vq actually converges strongly to Vq°° in L ( 



'2tTU>d\ 



that 



Using a cut-off technique as in the proof of Lemma ^thereafter, we deduce from (5.5) 



,oo vT^oo 



AVq°^ • Vq' 



f ■ V(Z°°. 



(5.11) 



The weak convergence of Vg to Vg°° implies that the right-hand side of (5.6) con- 



verges to the right-hand side of (|5.11|). Consequently, 



AVq"" • Vg^ 



In 



(5.12) 



and, denoting by As the symmetric part of A, (5.12) is equivalent to 



Jin 



AsVq'^-Vq'^. (5.13) 

1/2 

As is of course a uniformly coercive tensor field, we can thus define its square root As . 

\Al/'Vq'^h2(^,y (5.14) 



It follows from (5.13) that 

W^s l/]vVg ||l2 



,1/2 



On the other hand, multiplying (5.10) by As , we obtain 



(5.15) 



Because of the uniform convexity of L^(M ), it is well known that (5.14) and (5.15) 
imply 



Al^Hi^Vq^ -^ Al/^Vq°° strongly in L'^{R'^). 



Multiplying (5.16) by A^ , we finally have 



1/2 



l/„Vg^^Vg~ strongly in L2(R'^) 



(5.16) 
(5.17) 

n 
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Lemma 7. Let A be a tensor field from M'' to M'^^'^ such that there exist A > and A > 
such that 

V^ e M"', a.e in X e M°', AjCp < A{x)C ■ ^ and \A{x)^\ < A|C. 

Consider u G L^y^(M ) solving 



- div (^Vn) = in M'^, 



(5.18) 



T/ien n is constant. 



Proof. We define a smooth cut-off function x ^ C°°(M ) sucli tliat % = 1 in tlie ball Bji, 
X = in M.\B2R and \\Vx\\l^(w^) < 2/i?. 



Multiplying the first line of (5.18) by x^ ^-nd integrating by parts, we obtain 



/ AVu ■ (Vu) X = - ^Vii • (Vx) u. 

Jr^ jRd 

Using the Cauchy-Schwarz inequahty, this yields 



Br 



|V^P<^||Vx||io 



1/2 



iVnl 



< 



2A 
RX 



B2r\Br 

, 1/2 



1/2 



IVnl 



B2r\Br 



B2r\Br 

1/2 



(5.19) 



ul 



B2r\Br 



Defining 



UR 



\B2r\Br\ Jb2r\Br 



u, 



it is clear that u — ur is also a solution to (5.18) so that the above computations are valid 



for u — Ur. Since V(m — ur) = Vn, we deduce from (5.19) that 



ivnr < 



Br 



2A 



1/2 



|Vu| 



B2r\Br 



1/2 



\U-Ur\ 



B2r\Br 



(5.20) 



We next apply the Poincare-Wirtinger inequality to n — ur on B2r\Br. There exists 
a constant C{R) which depends only on R such that 



/ \u-ur\'^ <C{R) I \Vi 

■JB2r\Br J B2r\Br 



An easy scaling argument shows that C{R) is equal to R times the Poincare-Wirtinger 
constant on B2\Bi, so that there exists a constant C such that 



/. 



\u-ur\^ < CR 



B2r\Br 



Wu\ 



B2r\Br 



We deduce from ( |5.20[ ) and ( |5.21[ ) that 

/ |Vnp<^- 



IVnl 



B2r\Br 



(5.21) 



(5.22) 



Since Vn G L'^(M°'), the left-hand side of (5.22) converges to f^^ |Vnp when R — t- oo, 



and the right-hand side of (5.22) converges to 0. Then Vn = and n is a constant. D 
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Lemma 8. For d > 3, consider a Ij -periodic tensor field A such that there exist A > 
and A > such that 

V^ G M'^, a.e in x G R'^, Aj^P < A{x)^ • ^ and \A(x)C\ < A|^. 

Assume that the unit cell Q contains an inclusion D, the boundary of which has regularity 
(ji,M Jqj- some < /u < 1, and such that dist{D,dQ) > 0. Assume also that Ap^r is Holder 
continuous in D and in Q\D. 

Let f be a function in L?'{Q). 

There exists a unique solution u G Lf^^{W^) to 

- div (AVu) = div (1q/) in R'^, 
Vn G L^{R'^), lim u{x) = 0. ^^'^^^ 

|a;|— >oo 

Defining also uq the unique solution to 

r - div (AVuo) = div (1q/) in O, 
[uoGH^iO), 

where O is a bounded domain of R'^ containing Q and such that dist{dO,Q) > 1, there 
exists a constant K which depends only on A, A, /i, d, f and the Holder exponents, and 
not on the domain, such that for \x\ > 1, it holds 

K K 

\u{x)\ < -r-f^iiY, \Vu{x)\ < —r^. 

Proof. Let Go be the Green kernel associated with A with homogeneous Dirichlet boundary 
conditions on dO, uniquely defined by 

r -div(^VGo(-,y)) = 5y inO, 

and G be the Green kernel associated with A on M'^, unique solution to 

J - div(AVG(-, y)) = 8y in R'^, 

\ G{;y) G Wl^liR'') n H\R\B{y, 1)). 

We deduce from arguments stated in [6l Lemma 4.2] and relying on |10| Theorem 3.3], and 
on [51 Lemma 16] when A is Holder continuous and [13| Theorem 1.9] when A is piecewise 
Holder continuous, that there exists a constant K depending only on A, A, ^, d and the 
Holder exponents, and not on the domain, such that 

V(x,y) G O, |Vj,Go(x,y)| < -. -^, |V,Vj,Go(x, y)| < -. r^, (5.25) 

\x — y\°- ^ \x — yy 

V(x,2/) G R\ \VyG(x,y)\ < ^^, \V,VyGix,y)\ < , ^ ,^ . (5.26) 

\x — y\"- ^ \x — y\"- 
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It is well known that uq solution to (5.24) can be represented as 



uo{x) = / Go{x,y)dW{lQf){y)dy. 
Jo 

It is also clear that the function u defined by 

u{x)= [ G{x,y)div{lQf){y)dy 

is a H^-^^OR.'^) function which satisfies 

-div (AVu) = div (1q/) 
in the sense of distributions. 



(5.27) 



(5.28) 



Integrating by parts in (5.27) and (5.28) for x ^ Q, we have 

uq{x) = / VyGo{x,y) ■ f{y)dy, u{x) = / VyG{x,y) ■ f{y)dy, 



(5.29) 



and then 



Vuoix) = I V,VyGoix,y) ■ f{y)dy, Vu{x) = I V,VyG{x,y) ■ f{y)dy. (5.30) 
JQ JQ 



Using estimates (5.25) and (5.26) in (5.29) and (5.30) respectively, we find that there 



exists a constant K depending only on A, A, /i, d, f and the Holder exponents, and not on 
the domain, such that for Ixl > 1, we have 



K K 

\uq{x)\ < and \Vuq{x)\ < —^ 



\u[x)\ < 



K 



x\ 
K 



and \Vu{x)\ < , ,^. 



(5.31) 
(5.32) 



The function u being in H^^^iW^), we deduce from (5.32) that Vtt G Lp'iW^). Conse- 
quently, u solves 



div(^Vu) = div(lQ/) in 



Vu E L 



"i na>d\ 



(5.33) 



We know from Lemma ^ that ( 5.33| ) has a solution unique up to an additive constant. 
It follows from (5.32) that u converges to zero at infinity, so that u = u unique solution to 



(5.23). 
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